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USE  OF  A  ZONE  MODEL  FOR  VALIDATION  OF  A  HORIZONTAL 
CEILING/FLOOR  VENT  ALGORITHM 


INTRODUCTION 

An  algorithm  has  been  developed  to  describe  the  exchange 
flow  through  a  horizontal  vent,  such  as  a  hatch  in  a  ship  [1]. 
Until  recently,  there  was  no  way  to  fully  characterize  the  flow 
through  a  horizontal  vent  (Figure  1) .  The  Bernoulli  equation 
has  historically  been  used  to  calculate  uni-directional  flow 
through  horizontal  vents.  The  Bernoulli  equation,  which  is 
valid  for  all  vertical  vent  (i.e.,  doors  and  windows)  scenarios, 
states  that  the  mass  flow  rate  through  the  vent  is  proportional 
to  the  square  root  of  the  difference  between  the  gas  pressure 
inside  the  compartment  and  the  surrounding  atmosphere.  Use  of 
this  equation  to  calculate  horizontal  vent  flow  leads  to  an 
anomaly  when  the  difference  between  the  interior  and  exterior 
pressure  of  the  compartment  reaches  zero.  The  Bernoulli 
equation  predicts  no  flow  which  is  clearly  impossible  when  the 
gas  above  the  vent  is  more  dense  than  the  gas  inside  the 
compartment.  Gravity  dictates  that  the  heavier  gas  change 
places  with  the  lighter  compartment  gas  and  hence  an  exchange 
flow  must  exist. 

The  validity  of  this  new  algorithm  was  tested  by 
incorporation  into  a  computer  model  written  to  simulate  a 
compartment  with  a  fire  inside,  which  is  ventilated  through  a 
hole  in  the  ceiling.  The  model  predictions  were  then  compared 
to  the  experimental  results  of  reduced  scale  fire  tests  [2]. 


MATHEMATICAL  MODEL  -  AN  OVERVIEW 

The  mathematical  model  was  written  to  simulate  a 
compartment  fire  ventilated  only  through  an  opening  in  the 
ceiling.  The  zone  modeling  approach,  as  opposed  to  field 
modeling,  was  used.  The  latter  describes  a  system  using  the 
basic  conservation  equations  of  mass,  momentum,  and  energy  in 
partial  differential  form.  This  method  provides  very  high 
resolution  and  detail  but  is  computationally  intensive  for  even 
the  simplest  fire  scenarios.  A  zone  model,  on  the  other  hand, 
can  reasonably  approximate  a  much  more  complicated  fire  scenario 
and  still  have  a  significantly  lower  computational  requirement. 
The  system  is  divided  into  zones  which  are  assumed  to  be  of 
uniform  temperature  and  composition.  Mass  and  energy  balances 
are  written  around  the  zones  to  describe  the  interactions 
between  them.  The  system  in  this  particular  case  was  divided 
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into  a  total  of  seven  zones:  the  liquid  fuel  surface,  lower 
flame  zone,  upper  flame  zone,  enclosure  walls,  glass  viewing 
plate,  ceiling  vent,  and  the  interior  gaseous  medium  located 
between  the  flame  boundaries  and  the  enclosure  walls  [3].  The 
system  is  shown  in  Figure  2 . 

Transient  phenomena  are  described  through  a  sequence  of  time 
steps.  During  each  time  step,  there  is  no  accumulation  of 
energy  within  the  lower  and  upper  flame  zones,  or  the  fuel 
surface.  The  change  in  mass  and  energy  of  the  interior  gaseous 
medium  is  assumed  to  be  constant.  Mass  and  energy  balances, 
along  with  the  equation  of  state,  provide  five  non-linear 
algebraic  equations  which  are  simultaneously  solved  for:  mass 
loss  rate,  lower  flame  entrainment,  upper  flame  entrainment, 
interior  gas  temperature,  and  net  mass  out  of  the  compartment. 
One  dimensional  unsteady  state  heat  conduction  with  convection 
boundary  conditions  describes  the  heat  transfer  through  the 
compartment  walls.  Numerical  methods  are  used  to  solve  the 
resulting  differential  equations  for  the  interior  and  exterior 
wall  temperatures.  The  equations  and  relevant  assumptions, 
which  describe  the  individual  zones  and  their  interactions  with 
each  other  at  any  given  time  step,  will  be  explained  in  more 
detail  in  the  following  sections. 


Fuel  Source  Zone 

The  liquid  fuel  is  contained  in  a  circular  fuel  pan  and  the 
fuel  surface  is  assumed  to  remain  level  with  the  rim  of  the  pan. 
Lip  exposure  is  known  to  affect  the  burning  rate  of  a  liquid  pan 
fire  [4].  The  temperature  of  the  fuel  surface  is  taken  to  be 
equal  to  the  normal  boiling  point  of  the  liquid  fuel  and  there 
is  no  conduction  through  the  fuel  or  fuel  pan  [5].  The  mass 
loss  rate  can  then  be  calculated  using  the  following  steady 
state  conservation  of  energy  equation  written  around  the  liquid 
fuel  surface: 


mvCllq(Tg-Tllq)  +  mvLHv  -  Aehe(.T1[  Te)  Qt,B 


The  amount  of  energy  needed  to  raise  the  temperature  of  the 
fuel  to  its  boiling  point,  and  then  to  vaporize  it,  is  equal  to 
the  convective  energy  supplied  to  the  fuel  surface  from  the 
flame  and  the  net  radiant  energy  from  the  surrounding  zones. 

The  specific  heat  of  the  liquid  fuel  is  assumed  to  remain 
constant. 
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Flame  Zones 


The  division  of  the  visible  flame  into  an  upper  zone  and  a 
lower  zone  is  adopted  from  the  work  of  Steward,  Fang,  Heskestad, 
Delichatsios  and  McCaffrey  [6-11].  There  is  fairly  good 
agreement  between  the  various  correlations  with  respect  to  flame 
height  and  temperature.  There  is,  however,  a  large  variation  in 
mass  entrainment  predictions  using  the  available  correlations 
[12].  Due  to  this  lack  of  consensus,  flame  height  and 
temperature  are  determined  using  McCaffrey’s  correlations  [11] 
and  an  "effective"  mass  entrainment  is  calculated  based  on  an 
energy  balance  around  the  flame  zones. 

The  lower  flame  zone  is  the  continuous  region  located 
directly  above  the  fuel  surface  and  is  characterized  by  a 
constant  centerline  temperature.  The  upper  flame  zone  is 
intermittent  in  nature.  It  is  the  pulsating  portion  of  the 
flame  which  raises  periodically  from  the  continuous  zone  at  a 
consistent  frequency.  The  centerline  temperature  in  this  upper 
flame  zone  decreases  with  height  above  the  fuel  surface  as  the 
entrained  air  cools  the  product  gases. 

It  is  assumed  that  most  of  the  combustion  is  occurring  in 
the  lower  flame  zone.  McCaffrey  found  that  the  energy  generated 
in  the  intermittent  regime  is  equal  to  the  depletion  of  energy 
through  radiation  [11].  Figures  3  through  5  clearly  show  that 
the  net  radiant  heat  loss  from  the  upper  flame  zone  is  a 
fraction  of  the  heat  generated  in  the  lower  flame  zone  due  to 
combustion  for  all  cases  of  interest. 

The  steady  state  conservation  of  energy  equation  as  applied 
to  the  lower  and  upper  flame  zones  are  as  follows: 

Lover  Flame  Zone 


mvCvTe  +  mg_lfcgTg  +  -  (mg> lf+mv)  CltTlt  +  Ajxe{Tlt-Te)  +  0Iilf 


The  enthalpy  of  the  vaporized  fuel  and  the  entrained  air, 
and  the  energy  from  the  combustion  of  the  fuel  must  equal  the 
enthalpy  of  the  lower  flame  zone  products,  the  energy  convected 
to  the  fuel  surface,  and  the  net  radiant  heat  loss  from  the 
lower  flame  zone. 

Upper  Flame  Zone 


(mglf+mv)  C]fT 


it*  it 


mg.  UfCgTg 


<® g.lf+mg.uf+mv)  CutTtop  *  Qr.uf 


3 


The  enthalpy  of  the  lower  flame  products  and  the  entrained 
air  into  the  upper  flame  zone  must  equal  the  upper  flame  zone 
products  and  the  net  radiant  heat  loss  from  the  upper  flame 
zone. 


The  gaseous  specific  heats  are  assumed  to  be  constant  and 
equal  for  the  air,  fuel,  and  combustion  products  leaving  the 
visible  flame. 

The  products  leave  the  respective  zones  at  a  temperature 
based  on  an  assumed  squared  Gaussian  form  for  the  radial 
variation  of  the  temperature  [10,11]. 

Interior  Gas  Medium 

As  the  temperature  of  the  gas  within  the  enclosure 
increases,  there  must  be  a  net  loss  of  mass  through  the  vent 
during  each  time  step.  This  quantity  can  be  determined  by 
assuming  that  the  gas  within  the  enclosure  is  a  perfect  gas  and 
applying  the  equation  of  state. 


Once  sufficient  mass  loss  has  occurred,  the  pressure 
difference  across  the  vent  becomes  zero.  An  instability  is 
created  due  to  the  heavier  gas  above  the  less  dense  enclosure 
gas.  An  exchange  flow  occurs  which  is  calculated  using  the 
algorithms  described  in  reference  1. 


V  -  0 . 1 0  [  2  gAj  ' 5  (  p  amb-  p  g)  /  (  P  ?+  P  ami,)  1  0  ‘ 5 


The  total  mass  out  of  the  enclosure  is  then  calculated  by 
summing  the  appropriate  mass  flow  rates. 


m 


out 


-  +  m, 


'out,  not 
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The  rate  of  energy  change  during  each  time  step  is  assumed 
to  be  constant  and  can  be  described  by  the  following  equation: 


Eg  ~  ^g.  uf+^v)  ^UfT tDp  +  m  inC 'ajniJT 'amb  ^g,  lf+^g,  uf+^out^  CgTg  i?r<< 


hgl^gi  (  Tg  Tamb)  -  huAw  (  Tg  TJ) 


The  rate  of  energy  change  is  equal  to  the  enthalpy  of  the  flame 
products  and  of  the  air  entering  the  enclosure  through  the  vent 
minus  the  enthalpy  of  the  entrained  air  and  of  the  gas  exiting 
the  enclosure  through  the  vent,  the  net  radiant  heat  loss  of  the 
enclosure  gas,  and  the  heat  convected  to  the  walls  and  the  glass 
viewing  plate. 

The  temperature  of  the  enclosure  gas  is  determined  by: 


where  the  total  energy  of  the  gas  is: 


Eg  ”  Egi0ld  + 


EgUt) 


and  the  mass  of  the  gas  in  the  chamber  is  given  by: 


mg  ~  mg,old  * 


(itiv  mout  net)  At 


The  composition  of  the  enclosure  gas  is  recalculated  during 
each  time  interval  by  accounting  for  the  production  and 
consumption  of  gas  molecules  as  well  as  the  molecules  entering 
and  exiting  through  the  vent  opening. 

y3.«wt  °2,io  °2 . 
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^  HjO,  old  ™  -''jO,  on :  +  ™H2o.  io  ^H2o.pzoduc*d 


mN2  “  mN2.old  WN2.out  +  mN2,iD 


The  carbon  dioxide  and  water  concentrations  are  used  to 
calculate  the  absorption-emission  properties  of  the  enclosure 
gas  (described  below) .  The  average  molecular  weight  of  the 
enclosure  gas  is  required  to  calculate  the  net  mass  out  of  the 
enclosure  during  each  time  interval.  Analysis  of  the  enclosure 
gas,  on  a  dry  gas  basis,  shows  very  little  carbon  monoxide. 

Mass  Fraction 

Oxygen  0.14  -  0.23 

Carbon  Dioxide  0  -  0.08 


Carbon  Monoxide 


0  -  0.002 


Nitrogen 


balance 


Th»  ''ombustior  T-eacti^ns  are,  therefore,  assumed  to  go  to 
completion,  leaving  the  following  chemical  species  for 
consideration:  nitrogen,  oxygen,  carbon  dioxide,  and  water. 


Enclosure  Walls 


The  enclosure  in  this  model  is  comprised  of  the  walls, 
floor,  and  ceiling.  Energy  is  transferred  to  the  interior 
enclosure  via  radiation  and  convection.  One  dimensional 
unsteady  state  conduction  with  convection  boundary  conditions 
describes  the  heat  flow  through  the  enclosure  and  to  the 
surrounding  atmosphere.  The  equations  are  solved  numerically 
using  an  explicit  forward  differencing  technique  [13].  The  time 
and  distance  increments  were  chosen  based  on  stability  and 
convergence  criteria  as  well  as  optimization  of  computation  time 
versus  solution  accuracy. 

RADIANT  HEAT  TRANSFER  EQUATIONS 

Analysis  of  the  net  radiant  heat  exchange  between  the  zones 
is  by  the  net  radiation  method  [14].  The  surface  radiosities 
(enclosure  and  fuel  surface)  are  first  determined  by: 

Ji  -  <l-ei> +  eieu- 
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These  are  then  used  to  calculate  the  net  radiant  heat  loss  for 
each  surface: 

£ri  -  (Jt-GjAi 

and  each  gas: 

0ig  -  E  Ai*g(ebg~JJ 

Specifically, 

Fuel  surface: 

Or. s  -  U.-GJA. 

where : 

Jg  “  ®  s^bs  *  d~6B)Gfl, 

and 

GS  -  *  eJfebJf  +  tgT  Jfxu^si(l-Ffo) 

*  t  ufFsi  *  eufebufFgix  lf  +  Fait  ift  ufe  0e*,Ff0x  g 

Interior  Wall  Surface: 

where: 

l/j,  “  6  +  (  1  “®  GW 
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and 


l*oebo*gll  ~Ffo 


T  uf  ,  e  lfeblf T  ?  ^  ^  ^fo> 


Lower  Flame  Zone: 

0T.lf  “  ('As  +  Alf+Ai'>  e  lteblf  ~  Ase  lf^a  ~  Alfe  lfi  gebg  ~  Alf€ltXgJw 


-  Ap  1{eu{ebut  -  AjejfXufige^  -  Ate  lfx  ufx  gJ’w(l-Ffo)  -  A^e  1{x  ufc  0eboF{0x  g 

Upper  Flame  Zone: 

Qx '  uf  "  (Ai+Au[)e  u[ebu[  -  Ase  ute  lfebl[  -  A^  ufJBx  ieFgi  -  Aufe  u[e 

“  Auf&  uf  ^  g^w  ^  -  ^ lo^  Auf^  ut^  o^bo^foX  g 
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Enclosure  Gas: 


£r,s 


-  ( Au+Alf ■ 


uA-uf+^c)  ®  gebg 


-  gjy  ~  e  g^lfG  lfeblf  e  g^  uf^u(eu 


-  e  gAlf'Z  -  e  gAufl:  lfX  uf^a  A  Fsi  Aoecebo 

nuf 


Auf& 
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X  u lfeblf 


uf 


Initial  calculation  of  the  surface  radiosities  requires 
knowledge  of  the  view  factors  between  the  various  surfaces  and 
gas  boundaries.  The  view  factor  for  the  fuel  surface  to  the 
interface  between  the  lower  and  upper  flame  zones  must  be 
calculated  along  with  the  view  factor  for  the  visible  flame  to 
the  vent  opening.  The  remaining  view  factors  are  either  zero  or 
unity,  or  can  be  derived  using  the  reciprocity  theorem  or  the 
summation  rule  for  an  enclosure. 

The  view  factor  configuration  for  the  interface  between  the 
two  visible  flame  zones  and  the  fuel  surface  is  represented  as 
two  parallel  circular  disks  with  centers  along  the  same  normal 
[15].  For  the  purposes  of  the  radiation  calculations,  the  two 
zones  of  the  visible  flame  are  depicted  as  two  vertically 
aligned  cylinders.  The  view  factor  for  the  circular  vent  to  the 
visible  flame  is,  therefore,  also  determined  using  the 
configuration  of  two  parallel  circular  disks  with  centers  along 
the  same  normal. 

The  view  factor  between  the  square  vent  and  visible  flame 
should  ideally  be  calculated  using  a  circle  and  a  rectangle 
configuration.  Since  there  is  no  analytical  solution  for  this 
particular  geometry  [16],  the  "corners"  of  the  square  were 
ignored  and  the  previously  mentioned  circle  to  circle  geometry 
was  assumed.  This  translates  into  an  error  in  the  view  factor 
value  of  less  than  four  percent  for  a  typical  model  simulation. 
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The  heights  of  the  cylinders  were  calculated  using 
McCaffrey's  correlation  [10,11]: 


zu{  -  0 . 2  (iriyAH^  0,4 


zlf  -  0.08  (/n^)0-4 


Values  for  the  surface  and  gas  emissivities  are  required  to 
perform  the  radiation  calculations.  The  walls,  fuel  surface, 
and  glass  viewing  plate  are  assumed  to  be  gray  and  opaque. 

Their  emissivities  are  constant  and  equal  to  0.9,  0.42,  and  0.9, 
respectively  [3,13].  The  gases  within  the  two  visible  flame 
zones  and  the  enclosure  are  gray  and  transmitting.  The 
emissivity  of  the  gases  within  the  visible  flame  zones  are 
determined  during  each  time  interval  using  the  mean  beam  length 
approximation : 


e 
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1  - 


e 


-aL. 


where 


t.  - 

\  Area  / 


and  "a"  is  the  absorption  coefficient. 

Leckner's  approximation  [17]  is  used  to  calculate  the  emissivity 
of  the  enclosure  gas.  The  emissivity  is  a  function  of  the 
partial  pressure  of  H20  and  C02  within  the  enclosure.  The  mean 
beam  length,  the  temperature,  and  pressure  of  the  gas  are  also 
used  to  calculate  the  emissivity  of  the  gas  during  each  time 
interval . 


EXPERIMENTAL  STUDY 

Reduced  scale  fire  tests  were  conducted  to  obtain  results 
which  were  then  compared  to  the  model's  predictions.  The  test 
apparatus  consisted  of  a  0.43  m  (17  in.)  cube  with  0.0254  m  (1 
in.)  Kaowool  walls  and  0.0381  m  (1.5  in.)  Marinite  floor  and 
ceiling.  There  were  several  ceiling  boards,  each  with  a 
different  size  circular  opening,  which  allowed  for  variation  of 
the  vent  size.  A  0.216  m  (8.5  in.)  square  pyrex  plate  was 
positioned  on  one  wall  for  viewing  the  compartment  interior. 
Ethanol  was  burned  in  a  0.11  m  (4.25  in.)  fuel  pan.  The  fuel 
supply  system  was  designed  to  maintain  the  fuel  surface  flush 
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with  the  top  of  the  pan  rim,  thereby  eliminating  an  exposed  rim 
which  is  known  to  affect  the  burning  rate.  An  experimental 
measurement  of  the  mass  loss  rate  was  determined  from  a  load 
cell  located  beneath  the  fuel  supply  reservoir.  Nineteen 
thermocouples  were  located  throughout  the  compartment.  Interior 
gas  temperatures  were  measured  at  various  levels  along  with  the 
interior  and  exterior  wall  temperatures.  An  interior  gas 
analysis  gave  oxygen,  carbon  monoxide,  and  carbon  dioxide 
concentrations.  Two  tests  were  run  at  each  of  the  following 
vent  diameters:  0.1016  m  (4  in.),  0.127  m  (5  in.),  and  0.1524  m 
(6  in.).  Experiments  were  also  conducted  in  the  "no  vent"  (top 
completely  closed)  and  "maximum  vent"  (top  completely  open) 
configuration.  The  0.1016  m  (4  in.)  opening  was  not  large 
enough  to  sustain  burning  so  the  fire  self-extinguished  after  a 
few  minutes.  Further  details  of  the  experimental  set-up  and 
procedure  are  described  in  reference  2 . 


RESULTS 

The  experimental  results  for  the  0.127  m  (5  in.)  vent 
experiments  (exp  14  and  19)  are  shown  along  with  the  model 
predictions  in  Figures  6  through  11.  The  0.1524  m  (6  in.) 
experiments  (exp  15  and  18)  and  open  top  experiments  (exp  16  and 
17)  are  depicted  in  Figures  12  through  17,  and  Figures  18 
through  23,  respectively.  The  Figures  are  time  plots  of:  mass 
loss  rate,  oxygen  concentration,  carbon  dioxide  concentration, 
interior  gas  temperature,  interior  wall  temperature,  and 
exterior  wall  temperature.  In  general,  there  was  very  good 
agreement  between  the  model  predictions  and  the  experimental 
results.  There  was  also  very  good  experimental  reproducibility 
as  evidenced  by  the  plots  of  mass  loss  rate  and  interior  gas 
temperature.  Note,  one  of  the  assumptions  made  in  the  model  i 
that  the  interior  is  well  mixed,  i.e.  the  interior  gas  is  one 
zone.  In  reality,  there  was  stratification  with  respect  to 
temperature.  There  was  also  stratification  with  respect  to 
oxygen  and  carbon  dioxide  concentration  as  is  evidenced  in  the 
Figures  7,  8,  13,  14,  19,  and  20.  The  gas  sampling  point  was 
near  the  top  of  the  compartment  in  experiments  14,  15,  and  16, 
and  near  the  bottom  of  the  compartment  in  experiments  17,  18, 
and  19.  The  wall  temperatures  were  measured  only  in  experiments 
17,  18,  and  19. 

The  excessive  amount  of  scatter  in  the  fuel  mass  loss  rate 
measurement  makes  it  difficult  to  distinguish  a  rate  change  with 
time  for  any  of  the  experiments,  although  the  model  predicted  a 
slight  increase  for  both  the  0.127  m  (5  in.)  and  0.1524  m  (6 
in.).  The  model  did,  however,  correctly  predict  that  higher 
burning  rates  corresponded  to  larger  vent  openings  and  the 
predicted  rates  were  within  the  scatter  of  the  measurement.  The 
0.1524  m  (6  in.)  experiments  showed  an  approximate  0.3  kg/hr 
burning  rate,  whereas  the  0.127  m  (5  in.)  demonstrated  about  0.2 
kg/hr.  The  mass  loss  rate  in  the  open  top  experiments  averaged 
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0.5  kg/hr  which,  according  to  literature  predictions  [18],  is 
equivalent  to  an  open  fire  burning  rate. 

The  average  oxygen  concentration  decreased  with  time  in 
both  the  0.1524  m  (6  in.)  and  the  0.127  m  (5  in.)  experiments  as 
predicted  by  the  model.  The  model  predicted  a  higher 
concentration  of  oxygen  for  the  0.1524  m  (6  in.)  configuration 
relative  to  the  0.127  m  (5  in.).  This  was  also  verified 
experimentally.  The  model  slightly  under  predicted  the  oxygen 
concentration  for  the  two  smaller  vent  openings.  The  model 
predictions  for  the  open  top  configuration,  however,  were  very 
close  to  experimental  results  which  remained  slightly  below 
normal  atmospheric  conditions  for  the  entire  test. 

As  the  oxygen  concentration  decreased,  the  carbon  dioxide 
concentration  increased,  in  both  the  0.127  m  (5  in.)  and  0.1524 
m  (6  in.)  experiments  and  corresponding  model  predictions.  The 
increases  and  their  relative  magnitudes  were  correctly  predicted 
by  the  model.  Again,  the  open  top  experiments  were  only 
slightly  elevated  above  normal  atmospheric  concentrations  for 
the  entire  test. 

The  average  interior  gas  temperatures  at  the  end  of  the 
test  period  in  the  0.1524  m  (6  in.)  and  0.127  m  (5  in.) 
experiments  were  significantly  higher  than  the  temperatures 
reached  in  the  open  top  experiments,  i.e.  600  K  versus  400  K. 

It  appears  that  the  smaller  vent  openings,  while  reducing  the 
burning  rates,  also  served  to  restrict  the  loss  of  energy 
through  the  hole.  In  addition,  more  ceiling  surface  area 
translated  into  more  radiation  feedback.  This  rate  of 
temperature  increase  was  also  seen  in  the  model  results, 
although  somewhat  underpredicted.  The  interior  wall 
temperatures  corresponded  to  the  interior  gas  temperatures  both 
experimentally  and  in  the  model  predictions.  The  exterior  wall 
temperatures  were  similar  in  all  three  experiments  as  predicted 
by  the  model. 

The  conditions  within  the  compartment  apparently  had 
opposing  influences  on  the  behavior  of  the  fire.  The  fires  in 
the  0.127  m  (5  in.)  and  0.1524  m  (6  in.)  vent  opening 
experiments  burned  significantly  slower  than  the  fires  in  the 
open  top  experiments  even  though  their  environment  was  much 
hotter.  The  offsetting  condition  was,  of  course,  the  vitiated 
enclosure  gas.  Another  example  is  given  by  the  result  that  the 
open  top  experiments  showed  a  mass  loss  rate  comparable  to  the 
burning  rate  of  a  fire  in  the  open.  The  enhancement  expected 
due  to  the  radiation  feedback  from  the  walls  was  not  apparent 
even  though  the  interior  gas  temperature  was  elevated  compared 
to  normal  atmospheric,  i.e.  400  K  versus  300  K.  Again,  the 
experimental  results  showed  that  the  oxygen  concentration  was 
slightly  lower  than  normal  atmospheric. 
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SUMMARY 


The  validity  of  the  horizontal  vent  algorithm  was  verified 
by  incorporation  into  a  mathematical  model  and  then 
demonstrating  that  the  model  predictions  were  comparable  to 
experimental  results.  The  intent  was  not  to  accurately  model  a 
fire  burning  in  a  vitiated  environment  -  a  complicated  and  still 
largely  undefined  phenomenon.  The  intent  was  rather  to  show 
that,  for  a  given  burning  rate,  the  algorithm  was  able  to 
predict  from  the  area  of  the  vent  opening  and  the  density  of  the 
interior  gas  (and  thereby,  composition  and  temperature  of  the 
gas),  a  mass  flow  which  correctly  accounted  for  the  fire's 
oxygen  consumption  and  energy  output.  The  horizontal  vent 
algorithm  accomplished  this  on  a  small  scale  basis.  To  further 
test  and  confirm  this  algorithm,  a  similar  comparison  should  be 
made  on  a  larger  scale. 
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NOMENCLATURE 


A  surface  area 

C  specific  heat 

eb  black  body  emissive  power 

g  acceleration  of  gravity 

G  irradiation  (total  radiation  incident  upon  a  surface  per 

unit  time  and  per  unit  area) 
h  convective  heat  transfer  coefficient 

aHc  heat  of  combustion 

lHv  latent  heat  of  vaporization 

J  radiosity  (total  radiation  which  leaves  a  surface  per 

unit  time  and  per  unit  area) 

m  mass 

m  mass  flow  rate 

M  molecular  weight 

Q  power 

R  universal  gas  constant 

P  pressure 

T  absolute  temperature 

At  time  interval 

V  volume 

z  height 

p  density 

e  emissivity 

x  transmissivity 
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View  Factors 


Fsi  surface  to  flame  zone  interface 

Flo  flame  to  vent  opening 

Subscripts 

amb  ambient 

g  interior  gas 

gl  glass 

i  interface  between  flame  zones 

If  lower  flame  zone 

Jig  liquid  fuel 

o  vent  opening 

old  previous  time  interval 

out  out 

r  radiant  energy 

s  fuel  surface 

top  top  of  visible  flame 

uf  upper  flame  zone 

v  fuel  vapor 

w  interior  wall  surface 
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Figure  1.  Vent  Configuration. 
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Figure  2.  Modeled  System. 
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Figure  3.  Intermittent  flame  zone  net  radiant 
heat  loss  vs.  heat  of  combustion  in  lower  flame 
zone.  (0.127  m  vent  opening) 
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Figure  4.  Intermittent  flame  zone  net  radiant 
heat  loss  vs.  heat  of  combustion  in  lower  flame 
zone.  (0.1524  m  vent  opening) 
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Figure  5.  Intermittent  flame  zone  net  radiant 
heat  loss  vs.  heat  of  combustion  in  lower  flame 
zone,  (open  top) 
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Figure  6.  Model  predictions  versus  results  from 
0.127  m  (5  in)  vent  opening  experiments. 
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Figure  7.  Model  predictions  versus  results  from 
0.127  m  (5  in)  vent  opening  experiments. 
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Figure  8.  Model  predictions  versus  results  from 
0.127  m  (5  in)  vent  opening  experiments. 
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Figure  9.  Model  predictions  versus  results  from 
0.127  m  (5  in)  vent  opening  experiments. 
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Figure  10.  Model  predictions  versus  results  from 
0.127  m  (5  in)  vent  opening  experiments. 
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Figure  11.  Model  predictions  versus  results  from 
0.127  m  (5  in)  vent  opening  experiments. 
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Figure  12.  Model  predictions  versus  results  from 
0.154  m  (6  in)  vent  opening  experiments. 
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Figure  13.  Model  predictions  versus  results  from 
0.154  m  (6  in)  vent  opening  experiments. 
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Figure  14.  Model  predictions  versus  results  from 
0.154  m  (6  in)  vent  opening  experiments. 
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Figure  15.  Model  predictions  versus  results  from 
0.154  m  (6  in)  vent  opening  experiments. 
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Figure  16.  Model  predictions  versus  results  from 
0.154  m  (6  in)  vent  opening  experiments. 
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Figure  19.  Model  predictions  versus  results  from 

open  top  experiments. 
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Figure  20.  Model  predictions  versus  results  from 

open  top  experiments. 
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Figure  21.  Model  predictions  versus  results  from 

open  top  experiments. 
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Figure  22.  Model  predictions  versus  results  from 

open  top  experiments. 
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Figure  23.  Model  predictions  versus  results  from 

open  top  experiments. 
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